Mon. Not. R. Astron. Soc. 000, [TV?? (2005) Printed 20 August 2010 (MN WT^i style file v2.2) 



AGN have Underweight Black Holes and Reach Eddington 



A. R. King 1 

1 Theoretical Astrophysics Group, University of Leicester, Leicester LEI 7RH 



o 

(N 
< 



O 

u 

43 
6 



> 

m 
cn 

O 

o 



13 



20 August 2010 



ABSTRACT 

Eddington outflows probably regulate the growth of supermassive black holes 
(SMBH) in AGN. I show that effect of the Rayleigh-Taylor instability on these outflows 
means that SMBH masses are likely to be a factor of a few below the M — a relation 
in AGN. This agrees with the suggestion by Batcheldor (2010) that the M — a relation 
defines an upper limit to the black hole mass. I further argue that observed AGN black 
holes must spend much of their lives accreting at the Eddington rate. This is already 
suggested by the low observed AGN fraction amongst all galaxies despite the need to 
grow to the masses required by the Soltan relation, and is reinforced by the suggested 
low SMBH masses. Most importantly, this is the simplest explanation of the recent 
discovery by Tombesi et al (2010a, b) of the widespread incidence of massive ultrafast 
X-ray outflows in a large sample of AGN. 

Key words: accretion: accretion discs - galaxies: formation - galaxies: active - black 
hole physics 



1 INTRODUCTION 

It is now well established that the centre of almost every 
galaxy contains a supermassive black hole (SMBH), and that 
in many cases the mass of this hole is closely connected with 
properties of the host galaxy bulge through the M — a and 
M — Mbuigo relations (Ferrarese & Merritt, 2000; Gebhardt 
et al. 2000; Haring & Rix 2004). The connection is physi- 
cally reasonable, since the black hole binding energy r\M(? 
considerably exceeds the bulge binding energy ~ Mb u i ge o" 2 
(King, 2003) (here 77 ~ 0.1 is the accretion efficiency and a 
the velocity dispersion of the bulge). The black hole com- 
municates its presence to the host by driving powerful out- 
flows when it is fed matter at super-Eddington rates. If these 
outflows are momentum-driven, i.e. communicate only their 
ram pressure to the surrounding interstellar gas, the M — a 
relation emerges naturally as specifying the black hole mass 
at which an Eddington outflow can drive a significant bub- 
ble into the bulge gas (King, 2003; 2005). At smaller masses, 
the black hole can only drive bubbles which recollapse, and 
evidently do not interrupt the gas supply to the black hole 
which ultimately powers its growth. 

This argument implicitly suggests that the M — a rela- 
tion is an upper limit to the black hole mass, rather than a 
tight relation. It is now clear (Batcheldor, 2010) that this is 
probably so, as observational selection makes it difficult to 
measure black hole masses below the relation (cf Section 3). 

The obvious question then is how far below this limit the 
majority of SMBH lie. New insight into this question comes 
from recent X-ray observations by Tombesi et al. (2010a, b) 



of fast (v ~ 0.1c) outflows in a large fraction of local AGN. 
I shall argue here that this means that most local AGN 
contain black holes lying below the M — a limit, and that 
most of these systems undergo super-Eddington episodes 
which switch off only for relatively short intervals. The hole 
masses are probably not very far below the M—a value. I will 
show that Eddington outflows at masses significantly below 
this are Rayleigh-Taylor unstable and therefore inefficient 
in suppressing accretion. 



2 BLACK HOLE GROWTH 

The argument by Soltan (1982) relates the mass density of 
black holes in the local universe to the total background ra- 
diation they produced while growing. It suggests that the 
average medium-to-large galaxy hosts a black hole of mass 
> 10 s Mq. It is reasonable to suppose that the growth phases 
of these SMBH are observable as AGN. But since the inci- 
dence of AGN among all galaxies is relatively low, this must 
mean that the holes spend much of their time growing at the 
maximum possible rate, i.e. that specified by the Eddington 
limit. In this case the M — a relation constrains their growth, 
and we should expect the black hole masses in AGN to lie 
below the relation, with a galaxy nucleus ceasing to be ac- 
tive once its SMBH reaches this mass. We are thus led to 
the conclusions that 

(a) AGN must radiate at close to the Eddington luminosity 
for much of their lives, and 



2 A. R. King 



(b) their black holes must be underweight, i.e below the M — a 
relation. 

Although these conclusions follow from accepting that 
observed AGN do represent the growth phases of SMBH, 
neither of them is widely adopted in practice. This amounts 
either to simply ignoring the argument, or to tacitly assum- 
ing that SMBH growth only happens in AGN which are 
somehow unobservable. I shall argue below that this tacit 
assumption cannot be correct. 



3 M - a AS AN UPPER MASS LIMIT 

A major reason for the reluctance in accepting conclusion (b) 
above has probably been the tightness of the observed M — a 
relation. However Batcheldor (2010) has recently pointed 
out that the SMBH masses populating the observed M — a 
relation are found by using the velocities of stars within the 
SMBH sphere of influence 



Rint = 



2GM 



(1) 



Since one must resolve _Ri n f in order to get an unambiguous 
dynamical mass estimate, measured values of M correlate 
with a. Batcheldor (2010) shows that it is unlikely that one 
can measure masses M significantly below the value given by 
the observed M — a relation. Accordingly we are at liberty to 
regard the observed relation as an upper limit on the SMBH 
mass for a given a. 



4 SMBH FEEDBACK AND THE M a 
RELATION 

We have inferred that SMBH grow at the Eddington rate 
during active phases. The resulting massive outflows are ob- 
served (see Section 5). They offer a natural way of com- 
municating some of the hole's binding energy to the host 
galaxy, and thus establishing the M — a relation as they 
shock against the host interstellar medium. Whatever the 
Eddington ratio rh — M / M^dd, the accretion luminosity 
is always ~ I/Edd, because radiation pressure expels excess 
accretion at each radius to ensure that the local Edding- 
ton limit is not exceeded. Shakura & Sunyaev (1973) show 
for example that in a disc geometry the luminosity limit is 
of order Lsdd(l + lnm), where L^dd is the standard expres- 
sion for the spherical Eddington limit. As the outflow results 
from radiation driving, there are two extreme possibilities (cf 
King, 2003; 2005; 2010). 

For rh ~ 1 most photons in the AGN radiation field 
scatter about once and transmit their total momentum 
^Edd/c to the outflow. Further, the outflow shocks against 
the host ISM are effectively cooled by Compton cooling from 
the AGN radiation field. Thus only the outflow ram pressure 
^Edd/c is communicated to the host: this is a momentum- 
driven outflow. The outflow sweeps up the host ISM in a 
shell, and solving the shell's equation of motion shows that 
this recollapses unless the SMBH mass has reached the crit- 
ical value 

M a (mom) = J^-a 4 (2) 
(where f g ~ 0.16 is the cosmic baryon fraction wrt dark 



matter). If the hole mass is > M a the shell can reach sig- 
nificant radii within the host and presumably suppress cen- 
tral accretion and SMBH growth. A simple way of deduc- 
ing this result is to equate the thrust I/Edd/c to the weight 
of the swept-up gas shell. For an isothermal density distri- 
bution the weight of a swept-up gas shell turns out to be 
W — 2f g a 4 /G, independent of radius R (the swept-up mass 
goes as 7?, while its gravity ~ GM{R) / R 2 goes as 1/R - see 
King, 2010). Hence setting Z/Edd — W gives an estimate of 
the mass. (The solution of the shell's equation of motion 
shows that the radiation thrust accelerates the outflow to a 
fixed terminal speed, justifying the approximation; cf King, 
2005.) 

At the opposite extreme, for rh » 1, multiple scatter- 
ing within the outflow means that almost all the total pho- 
ton energy is given to the outflow. The outflow shock against 
the ISM is now not effectively cooled, and the total energy 
rate LEdd now does PdV work against the weight of the 
swept-up interstellar gas. This is an energy-driven outflow, 
as originally considered by Silk & Rees (1998). Equating the 
energy injection rate I/Edd to the rate of working Wa as the 
shell begins to move to large radii we find 



M a (en) ~ -4^<r 5 



(3) 



jcG 2 c 

The relations (J2j [3]) have no free parameter. Each is po- 
tentially self-consistent: the accretion rate is likely to be 
close to Eddington for the high-mass momentum-driven 
sequence (i.e. rh ~ 1), and significantly super-Eddington 
for the low-mass energy-driven sequence (i.e. rh >> 1) (cf 
King, 2010). The momentum-driven value M a (mom) is very 
close to the observed relation 



M ~ 2 x lO 8 M o-2oo- 



(4) 



(Ferrarese & Merritt, 2000; Gebhardt et al, 2000). 

In contrast, the relation @ is a factor ~ a/c ~ 10 -3 
below the observed one, and would imply unobservably small 
hole masses. It is important to understand why this low- 
mass sequence is apparently disfavoured, as if not, it would 
constitute a bottleneck to SMBH growth, and make it hard 
to understand why SMBH apparently do reach the higher- 
mass sequence @. 

The reason why the energy-driven sequence does not 
seem to occur in nature appears to be that the very strong 
density contrast between the shocked wind and the host in- 
terstellar medium in the energy-driven case causes the bub- 
ble blown by the outflow to break up through the Rayleigh- 
Taylor instability. By contrast, in a momentum-driven out- 
flow the postshock density is far higher, making the cor- 
responding bubble Rayleigh-Taylor stable for masses near 
M a , and potentially able to cut off accretion. In this latter 
case, the speed v of the preshock outflow is given by the 
condition 



Mv = ^ 
( 

as 

n 

V — — c 
m 



: rjM E ddC 



The continuity equation gives the preshock density as 



M 



_ rh M E dd 
4rrbR 2 v 4irR 2 br]c 



(5) 



(6) 



(7) 



AGN Black Holes: Underweight and Eddington 3 



at radius R, where b is the fractional solid angle of the out- 
flow (we shall see later that b is of order unity). As the 
outflow shock is efficiently cooled, the density is strongly 
increased (by a factor ~ { v / a ) 2 ) there. Thus the outflow 
density in contact with the host ISM is 



77777 M E ddC 



4-rvR 2 ba 2 

The ISM density is roughly the isothermal value 



PISM 



f 9 v 2 



2tyR 2 G 
so combining we find 



PISM 
Pm 



b M a (mom) 



2m 4 



M 



(8) 



(9) 



(10) 



Thus the bubble is Rayleigh-Taylor stable, and can propa- 
gate outwards and suppress accretion provided that 



(11) 



Since b ~ rh ~ 1, we see that central accretion and SMBH 
growth is likely to be suppressed for M ~ M CT (mom). (Note 
that accretion of gas within the black hole sphere of influ- 
ence can continue for a time after the outflow bubble begins 
to propagate to significant radii.) This suggests that SMBH 
masses in AGN are likely to be only a factor of a few below 
the M — <r limit ([2]), as they can grow only when the exter- 
nally imposed accretion rate is not far above the Eddington 
value and so potentially spend longer at such masses. 

For an energy-driven outflow, the preshock outflow ve- 
locity follows from 



T^Mv 2 = Z/Edd = r/M E ddC 2 



1/2 



(12) 



(13) 



With this change, and compression only by a factor 4 in the 
outflow shock (since cooling is negligible), the SMBH growth 
suppression condition corresponding to (|ll|l becomes 



M > 



2m 3 /- 



M CT (en) ~ 10 9 



7,3/2 



M CT (en) 



(14) 



which is clearly impossible to satisfy for realistic parameters. 
Thus an energy-driven bubble is never Rayleigh-Taylor sta- 
ble and cannot halt SMBH growth. We conclude that SMBH 
growth is not halted at the energy-driven sequence , and 
given an adequate mass supply can continue all the way to 
the momentum-driven limit ([2]). 



5 OBSERVED OUTFLOWS 

We have seen that SMBH feedback occurs through massive 
outflows carrying the Eddington momentum. These outflows 
must have speeds v ~ 0.1c (cf eqn[6}. Moreover since the 
momentum outflow rate is specified by the quantity pR 2 v 
(cf eqn[7|, the ionization parameter 



NR 2 



hLEddrrip 

p~R2 



(15) 



is also specified (here Li = hL^dd is the ionizing luminos- 
ity of the AGN, with U given by the spectral shape and 
ionization threshold). King (2010) shows that this condition 
requires that the outflows should have ionization parameters 

i = 3 x loSg.x^m" 2 , (16) 

where h = Zi/10~ 2 and 770.1 = 77/0. 1 , and thus show lines 
in the X-ray region of the spectrum. This is indeed what is 
observed (Pounds et al., 2003a, b; O'Brien et al, 2005). This 
reasoning shows that any outflow with velocities ~ 0.1c seen 
in X-ray lines carries mass and momentum rates of order 

A^Edd, ^Edd/c. 

The papers by Tombesi et al. (2010a,b) show that out- 
flows with these properties are extremely common. They 
are detected in a significant fraction (> 35%) of a sample 
of more than 50 local AGN. This must mean first that the 
solid angle factor b cannot be small (Tombesi et al. deduce 
b > 0.6), and second, that a large number of local AGN have 
undergone or continue to undergo episodes of Eddington ac- 
cretion. 

We get more information from the hydrogen column 
densities, which are observed to lie in the range 10 22 cm -2 < 
N H < 10 24 cm -2 . Using eqn. (0 we find 



N H = 



^dR- 



R in 



m 2 M Ed d _ 3 x 10 37 7h 3 M 8 cm -2 (17) 
AirbncmpRin brfc A Riu 



where _Ri n is the inner radius of the flow. For a continuing 
outflow this would be of order 100 Schwarzschild radii (corre- 
sponding to the escape speed ~ 0.1c). It is difficult to detect 
the corresponding very high column densities ~ 10 24 cm -2 
as the gas is likely to be fully ionized. The highest observed 
columns presumably correspond to continuing steady out- 
flows which have recombined at some distance from the black 
hole. Lower columns may reveal cases where the outflow is 
episodic, and last stopped at a time t g before the observa- 
tion. In this case R ln takes a larger value ~ vt g. Using JT7} 
we get 



t oS = 0.2 



m A M s 



yr, 



(18) 



where M 8 = M/lO 8 M ,iV 2 3 = AV(10 23 cm -2 ). Thus even 
the lowest observed columns (~ 10 22 cm -2 ) correspond to 
outflows which switched off only 2 years ago. This suggests 
that outflows are even more common than one might expect 
from the simplest interpretation of the results of Tombesi et 
al (2010a, b). 



5.1 Eddington or not? 

Tombesi et al. (2010b) consider three sources in detail (3C 
111, 3C 120, 3C 390.3) and suggest that their luminosities 
are below L E dd- However this procedure uses black hole 
masses based on assumptions that either the M — a rela- 
tion, or the SMBH-bulge mass relation, hold for all SMBH, 
at least in a statistical sense. (The M — a relation is used to 
calibrate the reverberation masses for 3C 120, 3C 390.3; Pe- 
terson et al., 2004, and the SMBH mass for 3C 111 uses the 
stellar bulge luminosity; Marchesini et al., 2004.) We have 
seen that both relations are likely to give overestimates of 
the SMBH mass, and hence overestimates of LEdd- Figure 
16 of Peterson et al., (2004) shows that a reduction by a fac- 
tor ~ 3 of SMBH masses estimated in this way would put 
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most of the sample of 35 reverberation-mapped AGN close 
to their Eddington luminosities. In addition to this, estimat- 
ing whether an AGN is near LEdd requires us to know not 
only its SMBH mass, but also its true bolometric luminosity 
Lboi, both to high accuracy. The latter problem is unlikely 
to be easier than the former. 



6 CONCLUSIONS 

This paper has argued that the black hole mass is a factor 
of a few below the M — a mass in active galaxies, and that 
a large fraction of AGN are fed mass at a super-Eddington 
rate, accreting just the Eddington value and expelling the 
excess. 

The first point follows from noting that SMBH growth 
towards the momentum-driven limit © is inevitable given 
a sufficient mass supply. In particular, energy-driven out- 
flows are Rayleigh-Taylor unstable, so the mass is not con- 
strained by the energy-driven limit Q. Growth only slows 
when momentum-driven outflows become Rayleigh-Taylor 
stable, i.e. when the black hole is a factor of a few below 
the M — a value. So SMBH masses in AGN are likely to 
be below, but fairly close to, this critical value. This agrees 
with the suggestion by Batcheldor (2010) that the M — a 
relation is an upper limit to SMBH masses. 

The idea that AGN regularly reach LEdd follows natu- 
rally from noting that SMBH have to grow rapidly to reach 
the masses specified by the Soltan relation. It is consistent 
with the first proposition above (low SMBH masses): Ed- 
dington ratios for observed AGN must be higher than pre- 
viously estimated if their black holes lie below the M — a 
relation rather than on it, as is sometimes assumed. 

The strongest evidence for Eddington accretion comes 
from the papers by Tombesi et al. (2010a, b), which show 
that a large fraction of nearby AGN have outflow with ve- 
locities ~ 0.1c and ionization parameters £ ~ 10 4 — 10 5 , as 
expected. At face value this suggests that a large fraction 
of local AGN are fed at super-Eddington rates, and I have 
argued that it is difficult to avoid this conclusion. 
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